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Abstract 

O 
O 

The thermoemf in bipolar semiconductors is calculated. It is shown that it is 
necessary to take into account the nonequilibrium distribution of electron and 



hole concentrations (Fermi quasilevels of the electrons and holes). We find 

that electron and hole electric conductivities of contacts of semiconductor 

samples with connecting wires make a substantial contribution to thermoemf. 
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INTRODUCTION 



In calculating thermo-emf and explaining its nature, it is common to consider at first 
the unipolar case.0 The physical transparency of this phenomenon and the clearness of 
calculation lead to the following paradoxical results: in the case of a bipolar medium the 
situation seems to be equally obvious, and so the same calculation scheme is used.0 The aim 
of this paper is to show that the situation changes in bipolar media in principle, from the 
point of view of both the physics of the proceeding processes and the methods of thermo-emf 
calculations. 

Let us begin with the unipolar situation. Dating back to the paper of Thomson! in 1856 
the theory of the origin of the thermo-emf in this case has found its most precise description 
in the publications of Seeger and Kaydanov and Nuramski.i It is necessary to connect the 
electric circuit (Fig. [l|) for the determining thermo-emf in the case of the absence of electric 
current (broken circuit). A semiconductor sample whose thickness is 2a (—a < x < a) 
contacts with a heater with temperature T\ on the surface x = —a and with a cooler with 
temperature T2 on the surface x = +a. The connecting wires are made of the same material 
(metal) and are hooked up to the terminals of a measuring compensating circuit which allows 
us to measure the difference of voltage ipb and <p c in the absence of an electric current. Leads 
of connecting wires have equal temperatures [for example, T* = (7\ + T 2 )/2] at points b 
and c. As follows from Ohm's law, for a closed electric circuit V = j(R + r) (V is the emf of 
the power source, r is source's resistance, R is the r esistance of the external load, and j is 
the density of the electric current). If R — > 00 (broken circuit), then the density of electric 
current j — > and V = jR = ip c — (^(,.1 

Since the electron's chemical potentials n n (x) are equal in the points b and c, then 



(f c -cp b = {p, b n - fl c n )/e, 



(1) 



where fi n (x) = /U n (x) — eip(x) is an electrochemical potential. At the same time, 




(2) 
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As is well known, the expression for the density of the electric current is of the form§ 

j„ = vntfi) - a^VT^j , (3) 

where cr n (x) is the electric conductivity and a n (x) is the thermoelectric power. In the absence 
of an electric current djx n {x)/dx = ea n (x)dT/dx. Let us emphasize that this correlation is 
correct all over the circuit, where cr n (x) ^ 0, e.g., outside the region [b,c], where a n (x) — > 
(this condition provides j n = everywhere in the circuit). Therefore the expression (§) can 
be rewritten in the form 

fin ~ fin — / eot n {x)—dx= I eaMdT+ / ea n dT + eaudT, 
Jc dx Jt* Jt 2 jt x 

where a n and au are the values of a n (x) in the semiconductor and in the connecting wires. 
Finally we have for the difference ip c — (fbi which coincides with the thermo-emf V of the 
broken circuit, 

V = y? c - Lp b = (a n - a M ){Ti - T 2 ). (4) 

To simplify the calculations we have assumed that T\ — T2 T* which does not restrict 
the generality. Since usually |«m| <S \&n\, then 

V = a n (T 1 -T 2 ). (5) 

The last expression coincides formally with /_* V[y — (/i n /e)] dx. 

Thus the general scheme of calculating the thermo-emf comes to the following for a 
unipolar broken circuit. The gradient of the electrochemical potential Vfi n (x) is obtained 
by setting expression (|j) equal to zero. Then integration of S/fi n (x) in the anticlockwise 
direction (together with multiplying by e~ x ) gives the value of the thermo-emf. It is impor- 
tant to emphasize that, as above, the extreme points of the circuit (b and c) would have the 
same temperatures. 

It follows from (|2|), in a unipolar medium, that even if the condition of quasineutral- 
ity is not fulfilled and nonequilibrium electrons (this notion will be defined more exactly 



below) arise, the gradient of concentration of these carriers n does not give a contribu- 
tion to the thermo-emf [the summand (9// n /9n)Vn disappears from the expression for V]. 
This argument confirms the following idea of Ioffe:@ the potential difference created by the 
bulk concentration's gradient is compensated by the contact potential difference ( "diffusion 
voltages") on the boundaries +a and —a. 

Let us note in conclusion that the obviousness of the above scheme is lost when thermo- 
electric current flows in a closed circuit [j n ^ 0, see Eq. (|13|)]. 



I. BIPOLAR SEMICONDUCTORS. TRADITIONAL APPROACH 



The system of equations for electrons and holes are analogous to 



cr n (X 



Jn w n 



Jp = Vp{.X) 



+v Mx) _ (x)VT 



-V^p - a p (x)VT 



(6) 



Here j p is the density of the hole electric current, a p (x) and a p (x) are the electric con- 
ductivity and thermoelectric power of the holes [a n (x) and a p (x) have opposite signs], 



is 



fip = ~£g ~ fin = flp + e(p 

the hole electrochemical potential,^ and e„ is the band gap. It is important to emph 



(7) 



that in expressions (j^j V/i = (<9/i/<9T)VT usually. The full current j is equal to 



j = jn + j P = 0, 



asize 



(8) 



if a bipolar semiconductor is represented as in Fig. [I]. 

It is easy to obtain V[fi n (x)/e] from (j|) taking into account (|6]) and knowing that V/2 F 
-V/J n [see (0)], 

fin(x)] I 



alx. 



[a n (x)a n (x) + <j p (x)a p (x)) VT, a(x) = a n (x) + a p (x). 



Then we find the thermo-emf V (it is assumed that a n , p ^> au and T1—T2 <C T* as above): 



V 







1 c 


e 



dx = -(cr n a n + o p a p )(T 1 -T 2 ). (9) 
cr 



This expression, especially the method of the calculation, causes serious objections. 
Really, the electron concentration and hole concentration should be lower in the heating 
lead in the stationary state due to thermodiffusion if bulk and surface recombin ations 
are absent. In contrast, these concentrations should increase in the cooling lead. On the 
one hand, this causes the appearance of appreciable diffusion currents! in expression @ 
[V/i n = (<9/i n / 'dn)Vn\. On the other hand, it leads to a violation of relation (|7|). Two Fermi 
quasilevels jX n and \x v arise instead of a single level of the electrochemical potential, and as 
a result |V/2 n | ^ |V/t p |. In this case the procedure described in the beginning of this sec- 
tion becomes incorrect, because the single common "gradient of electrochemical potential" 
of electrons and holes is absent. Moreover, if bulk and surface recombinations are absent 
then both partial currents j n and j p should be equal to zero, not only the full current j 
(in + jp = 0) . As a result we have two equations (j n = 0, j p = 0) for both thermoelectric 
fields V(/i n /e) and V(/i p /e) instead of one equation (??). One more problem arises when 
bulk and surface recombinations take place: the correct determination of electron and hole 
equilibrium concentrations.0 

Finally, the question remains how to obtain the thermo-emf in this case, and which 
physical phenomena determine its value. The answer to the first of these questions has been 
given,! where the general scheme was proposed for calculation of an emf of any nature. It 
follows from this paper thatS 

v=f( ^ - dl-l (™± + ^dl, (10) 

J \ea dx ea dx J J V cr a J dx 

where integration is carried out clockwise. Let us note that expression fllQ|) is always correct 
(for a broken circuit just as in the case of a flowing thermoelectric current). The second 
item in (PH) coinci des with the expression (|9]). The first item in ([H]) vanishes (it is equal 
to zero identically in the unipolar case) and expression ( |T0D turns into @ if electrons and 
holes have a single, common electrochemical potential level. Correlation (0) does not hold 



if electrons and holes flow from the hot lead to the cool one and electron and hole Fermi 
quasilevels appear. In this situation the first item in (|T0| ) differs from zero and the gradients 
of concentrations and corresponding diffusion currents contribute to the thermo-emf. 



II. THERMOELECTRIC PHENOMENA IN BIPOLAR MEDIUMS 



Let us go on to the description of an approach to exploration of the thermo-emf, which 
does not involve either the contradictions or the incorrectness pointed out above. Note that 
some aspects of this approach have previously been expounded.BIl 

Let us restrict ourselves to the first case 



AT 



(ii) 



for simplicity. Here T% is the heater temperature at the point x = —a, T 2 is the condenser 
temperature at the point x = +a, T* = (T + T 2 )/2, and AT = Ti — T 2 . 

The condition when the temperature field of the quasiparticles (electrons, holes, and 
phonons) is common and is a linear function of coordinates has been obtained earlier.^ Let 
us assume that 



AT < T*. 



(12) 



Let an arbitrary semiconductor be defined by the function /x°(T*) which is founded from 
the condition of electroneutrality.B Then becomes a function of the coordinate /i°(T(a;)) 
in the temperature field ( [□]) , and 



M °(x)= M °(T*) + <5 M °(x), 5nl{x) 



X 



-AT- . 

dT* 2a' 



6 £ « /i°(T- 



(13) 



Function (|I3"D gives, uniquely, the concentration distribution in the sample: 

n (x) = n (T*) + 5n (x), 

~u°(T*) 

n (T ) = 7„(T )exp 
Sn (x) = n (T 



T* 



1\ ^L± + S A 

2 T* 2a T* 



(14) 



Here j n (T) oc T 3 / 2 is the density of states at the bottom of the conduction band. 
If we introduce (here it was assumed that the energy gap e g is independent of tempera- 
ture) 



Sr. 



AST 



then we can write 



analogously to ([13]). Then the hole concentration is 



(15) 



Po(x) =Po(T*) +8p {x), 



Po(T*) = 7 P (T*)exp 
8p (x) =p (T* 



3\ AT x 5fi° v 
1 + ' 



(16) 



T* 2 J T* 2a T* 

where 7 P (T) oc T 3 / 2 is the density of states at the top of the valence band. 

Let us note that condition (|12"D is not sufficient for the correctness of formulas (|l4f) and 
(|16|) in contrast to (|T3| ) and (|i5|) . The additional conditions 



Wn,p\ « T* 



AT 



max[|/i° |, 



are necessary. 

It is important to emphasize that n (x) and Po( a; ) are n °t "equilibrium" concentrations 
(see the Appendix) . Inverted commas are used here because it is impossible to use the term 
"equilibrium," strictly speaking, in the presence of a temperature field (|IID . 

The situation becomes nonequilibrated in the authentic sense (see the Appendix) when 
the gradient of electrochemical potential becomes nonzero because of taking into account 
the terms a n VT and a p VT. Let us examine this situation, assuming that bulk and surface 
recombinations are absent, and under the condition of a broken circuit (as in Sec. |). 

In this case the stationary distributions of concentrations 



7 



n(x) = rii(x) + 8n, p(x) = pi(x) + Sp 

and electric potential 

<p(x) = S(pi(x) + Sip 



are described by the system of equations j n = 0, j p = and 

= ^e(Sn - Sp) 

ax 2 

[expressions for ni(x), pi(x), and 8xpi{x) are contained in (|A2|) and (|A7|)1. 

Let us note [see @] that the quantities \p$ {ftp) are the "equilibrium" electron (hole) 
chemical potentials, see (|A2|) and (|A7 )1 



(-0 a (i) a 

ftn ftn vftni ftp ftp Sftp, 



5n = T>> 5ft n , 5p = ^ J 5 ft p , 



which are contained in the expressions for j n and j p are not already connected by expression 
(|H)|), e.g., they are compatible with Fermi quasilevels. As for quantities a n and a p , depending 
on ft n and ft p respectively^ they are determined by the quantities ft n (T*) and ft p (T*), since 
we carry out all calculations up to the members of order (AT/T*). Thus Sn Q , Spo, Sn\, 
5pi, and 5(fi (5ft n , 5 ft®, 5ft£\ Sftp, and 5pi) are not incorporated into the system of 
equations for finding Sn, 5p, and Sip (or 5ft n , 5ft p , and Sip). As a result, the determining 
of unknown quantities needs no calculations presented by formulas (??)-flT6|) and ( |A1| )- 



). We shall recall, however, that these calculations are necessary if we have to take into 
account recombination. Conditions ( |A5| ) and ip\ x=Q = are use d f° r determining integration 
constants. 
Then 

Sip = A x {e Xx - 1) + A 2 {e- Xx - 1) + qAT^, 
5 ft n = +eSip + e a n T(x) + e C\, 
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5/ip = —eSip — e a p T(x) + e C 2 . 



(17) 



Here 



a 



a n n (T*) + a pPo (T* 



n (T*) + p (T*) 

Constants A\ j2 and Cx t2 [see (|TT|)1 are connected by correlations 

no^d-poC^Ca 



n (T*) +p (^* 



A + A 2 - aT* 



£,pC 2 - inC\ = («p£ p + Ot n £ n )T 2 + (£„ + £p)<Pc 



[for definitions of £ n and £ p , see (|A4|) ] . 

For determining all the constants presented in ([17]) , it is necessary to give conditions 
5ip\ x=±a and 5fi n (or 8/i p ) at x = 0, for example.! 

Let [see (TAg)] 

Mx=±a = Pc,6, 5/Xn(0) = 5//p(0) = 0. 



Then expressions (|T7|) turn into 



5 ip 



AT 



x 

a 

a 



a + 



ctnjn + a p £ p \ sinh Xx 
€n + L J sinh Aa 



S^n = +e<5y? — ea n AT 



5 Hp = —e5ip + e a p AT 



x 

2a' 



(18) 



As follows from formulas ( |17|) and (|T8|), 5// n 7^ — <5/V e.g., two Fermi quasilevels really 
appear. 

If we use the condition of quasineutrality Aa ^> 1, then 



5(p = aAT 



x 

2a' 



5/x„ = e 7^- (a p - a n )AT—, 

n (T*) + po{T*) 2a 
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n (T*) x 
n {T*)+p (T *) 2a 

n (r> (T*) ATx 

on = op = e — 7-— — a B — a n )— . (19) 

y n (T*) +p {T*) K T* 2a K ' 

Thus there are two Fermi quasilevels even in the quasineutrality approximation and the 
condition of quasineutrality reduces to the equality Sn = 5p [compare with ([Sj])]. 

When the semiconductor is unipolar [for example, n Q (T*) 3> p (T*)] 5n = 5p <C 
Po(T*) <C n (T*) and <5« n = 0, e.g., the redistribution of concentration and Fermi level 
change do not take place. This is in accordance with the results of the Introduction. 

To conclude, let us note that expressions flT9| ) could be derived from the equations j„ = 
and j p = only without using the Poisson equation, assuming at once that the relation 
5n = dp takes p lace when Xa ^> 1. 



III. THERMO-EMF OF BIPOLAR SEMICONDUCTOR 

As was noted above, the thermo-emf is described by the expression (|10D in a bipolar 
medium, and it is important to emphasize that this expression does not contain the electric 
potential <p(x). The latter is quite natural if we wish to use the correct determination of 
the thermoemf, which is formed by forces of nonelectric origin, but not to use the artificial 
scheme which was presented in the first section. 

Let us assume that a semiconductor sample is placed in the interval —a < x < +a. It 
is connected with an instrument by metal wires with chemical potential um which does not 
depend on temperature. Let the thermoelectric power a M be equal to zero (o^f a n ,a p ). 
We assume that u n (T*) = fi M for simplicity. 

As was noted, the first integral in expression (|10D tends to zero when electrons and holes 
have a single level of chemical potential (5{i n = —5{i p ). So we have: 
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Taking into account that 8fi n and S/j, p oc AT [and so a n and a p depend on hq{T*) and 
Po(T*)], we get 

1 



V 



ea(T*) 



cx n (T*) r a d(6fi n )-a p (T*) f +a d(6fi p 

J— a J— a 



[5fi n (-a) - Sn n (+a)} - [Sfi p (-a) - 5/j p (+a)] 
1 Wrj^ + a^fjaJ^-r!); 



cr(T* / 

e = a + e P - (20) 

The second and third terms in the expression ([ZD) correspond to the contributions of 
the jumps of \x n and \i p on the surfaces x = ±a to the first integral (]T0|) . The analogous 
contribution of the second integral is absent, since in our problem temperature is continuous 
at x = ±a [compare with (|A5|) 1. Two additional terms [(Cn«„/0 + (£p a p/0]A^~± appear in 
expressions (|A5|) and (^) if the thermoconductivity of planes x = ±a is a finite quantity, and 
thus the temperature has discontinuities AT±.i The index "s" serves to show the relation 
to the planes x = ±a. 

We find from expressions (|20|) and (p!7|): 

V = ^ a "^ pQp AT. (21) 

For example, the Fermi quasilevel contribution to the thermo-emf compensates completely 
the conventional thermo-emf expression (Q). The nonzero thermo-emf is caused only by 
Fermi quasilevel breaks in the contacts of the semiconductor sample with connecting wires. 
Let us emphasize once again that the last assertion is true with any form of boundary 
conditions on planes x = ±a. 



CONCLUSION 



Thus the thermo-emf V is determined by thermoelectric powers a n , p and surface charac- 
teristics Comparing formula ( |2"TD with formula (Q) we see that taking into account Fermi 



quasilevels in the thermo-emf changes its value substantially. So if a n a n (T*) + a p a p (T*) = 

11 



then expression (|9|) becomes zero but the value V obtained from (|21| ) does not equal zero. 
In the general case expressions (|]) and (|2T|) may have different signs. 

Returning to the unipolar case (£„ 3> £ p ) we come to the expression (|5|) as was noted 
above. Let us notice only that the condition £ p 3> £ n does not follow from the condition 
Po(T*) ^> n (T*) in a hole semiconductor. The problem is that the ratio of the electron 
mobility to the hole mobility can become very large in some semiconductors [for example, 
it is more than 80 in InSb (Ref. g)], and the transition from expression ( fZlf) to @ does not 
always take place. 

It is clear that use of expression (|2l"D can essentially change the calculated value of the 
efficiency of a thermoconverter. We would like to finish with one comment. It is impossible to 
say anything about the order of magnitude of V in typical experimental situations, because 
nobody knows the value of £. 
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APPENDIX A: 

The chemical potential n° n {T(x)) [see ([13|) 1 corresponding to concentrations n (x) and 
Po{x) is heterogeneous in space. So diffusion currents will arise in the process of the estab- 
lishment of "thermodynamic equilibrium", which leads to redistribution of concentrations, 
the appearance of bulk charge layers, and an internal thermoelectric field characterized by 
the electric potential S(fi (x). Here the situation is analogous to the process of establish- 
ing thermodynamic equilibrium in heterogeneously doped semiconductors.0 The following 
distributions of the electric potential S(fi(x), concentrations, and chemical potentials: 

m(x) = n (x) + 8n x {x), p x (x) = p (x) + 5px(x), 
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fl W (z) = /i° (x) + 8^ (x) , //« (x) = //J (x) + fyW (x) , 

correspond to "equilibrium" when the electrochemical potential is constant: 

fJ>n( x ) ~ edipi(x) = —e g — ^p\x) — ed(pi(x) = const. (Al) 

The functions 8n\{x) and 5pi(x) are connected with 5fi$ and 5fj,f' [see ([b|) and (|16|)] 
naturally by the formulas 

Sn x = n (T*)^, 6 Pl = PoCH^b 8$> = (A2) 

Thus there are two unknown independent functions. It is necessary to use the Poisson 
equation to determine them, besides Eq. (|A1| ). 

d ^M- = A^e{8n 1 -8p l ). (A3) 

It is easy to formulate boundary conditions for Eqs. ([Al]) and ( |A3| ) if we introduce 
electron and hole electric conductances per unit area of the contacts of the semiconductor 
sample with connecting wires (£ n and £ p , respectively). If the thickness of the junction 5 
is negligible compared to the thickness 2a of the bulk semiconductor and <r* (cr*) is the 
junction conductivity of electrons (holes), we have (cr* and are supposed to depend on 

and ^=&„f' (A4) 

We choose the origin of <f(x) at the point x = [<^(x)| a , =0 = 0]. 

Then the condition of continuity of electric current in contacts in the broken circuit [con- 
sidering that a n = a p = in the connecting wires and hm does not depend on temperature 
(metal)] is reduced toi 



e e 



= (fn + Q<Pcb- (A5) 

x=±a 



Condition (K5I) in the "thermodynamic equilibrium state" turns into 
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x=±a 



0. 



since <p c = ip b = and S/i p = —Sfj, n . 

As a result, for 8y u 5^ (5/i^ = -8y$) and 8m [n (T*) 5 Pl = -p (T*) 8m] we find 



1 



8<p x = -5/ii 0) (x) + 2C sinh Ax, 

e 

8/1^ =2Ce sinh Ax, 
Sni = 2Ce^ sinhA*. 



(A6) 



Here 



A- 1 



rjiif 



1/2 



[4vre 2 [n (T*) +p (T*)] J 
is the Debye radius. 

It is necessary to formulate boundary conditions either on 8(fi or on 8jJ$ for determining 



the constant C. It is natural to assume that Sipi 



x=±a 



0. Then 



C 



dv n (T*) AT 1 



dT* 2e 2 sinhAa 

If the condition of quasineutrality holds, e.g., Xa 1, as usual, takes place, then 

<fr°(T*) AT x 
^ 1 = -^T^ — 2^' 

5m = = = ^(i) = o. (A7) 

At first sight it seems that Eq. (|A3| ) implies the condition 8m = 8p\ only. But we 
see from ( |A7| ) that the condition of quasineutrality reduces to the stronger requirement 
8n\ = 8px = in a bipolar medium as in a unipolar semiconductor during the process of 
establishing "equilibrium." 

It is clear that functions m{x) and p\(x) have to be named the "equilibrium" concentra- 
tions of the carriers because 



[ ^n\ x ) ~ e8(p x (x) =0, $\ x ) +eSip 1 (x) 



dx 



0. 



These functions must participate exactly as "equilibrium" concentrations in the expressions 
for the bulk and surface recombinations. 
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FIGURES 



FIG. 1. The electric circuit for measuring of thermo-emf. 
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